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Abstract 



A finite temperature many-particle theory of condensed matter systems 

is formulated using the functional Schrodinger picture. Using the interacting 

electron gas as a model system, we solve the equation of motion for the density 

matrix variationally with a Gaussian type trial density matrix. We show that 

the present formalism yields the finite temperature Hartree-Fock results both 

for the para- and ferromagnetic states in a simple and convenient fashion. 

Implications of the present results and future prospects are also discussed. 
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I. INTRODUCTION 



The many particle systems at finite temperatures are generally studied using the temper- 
ature Green's function through either diagrammatic or functional integral approach 
This Green's function method utilizes the Heisenberg picture and enables us to evaluate 
the equilibrium and nonequilibrium thermodynamic properties of systems such as the free 
energy and the linear responses when perturbative expansions are available. On the other 
hand, the Schrodinger picture, which is a more familiar concept in quantum mechanics has 
not been much used in the many-body theory, although the usefulness of this picture for 
the study of quantum structures of field theories was recognized early by several authors 
f§-§§. Only recently, it has been shown that a many-particle theory of condensed matter 
systems can be successfully formulated using the functional Schrodinger picture(FSP) 0. 
However, the theory was limited to the zero-temperature and, thus, requires extension to 
finite temperature to be truly useful for many-particle studies. 

In this paper, we present a finite temperature formulation of many-particle theory based 
on the functional Schrodinger picture. We apply the theory on the interacting electron gas 
system and show that a Gaussian approximation of the theory produces the finite tem- 
perature Hartree-Fock results both for the para- and ferromagnetic phases in a simple and 
convenient fashion. We also show that the self-consistent equation for the ferromagnetic 
splitting yields the Stoner criterion for the ferromagnetism. 



II. NON-INTERACTING ELECTRON GAS SYSTEM 

A finite temperature quantum field theory using the functional Schrodinger picture has 
been recently formulated MM by generalizing the Gaussian approximation of the variational 



calculation in quantum mechanics [10]. The method is based on the Feynman's observation 



that the unnormalized density matrix 

P = e-? H (1) 



2 



satisfies the equation 



dp rr 

dj3 =Hp 



(2) 



where (3 is the inverse temperature and H the Hamiltonian of the system [[□]] . 

First, we study the non- interacting particle system which has the grand Hamiltonian 
given by 



K = H — fxN 



x 



2m 



1pa(x) — PN 



(3) 



where N represents the number operator and /i the chemical potential. In the FSP, the 
fermion field operator ip and its conjugate momentum iifi are expressed as functional op- 



erators. Following the Floreanini-Jackiw prescription |12[], the fermionic field operators are 
expressed as 



ip(x) 



V2 



u(x) + 



u(x)^ + 



5 



5u(x) 



(4) 



to satisfy the equal-time anticommutation relation {ipi(x, t), ipj(x', £)} — 5ij5(x — x'), where 
u and w are anticommuting Grassmann variables. Using the above expressions, we rewrite 
the grand Hamiltonian as follows, 

5 \/ 5 



K = \Y.UbU a + 



A,B 



8u f 



u B 



(5) 



where ^ab = — vWfflaP ~ P^ab and A = (x, a), B = (y,f3). We choose a Gaussian 



trial density matrix 



p — e u\F AB u B +J A u A +u\L A +C 



(6) 



to satisfy the differential equation, 



d P is 



(7) 



Substituting Eq.(6) into Eq.(7), we obtain the following identity equation, 
,dF dJ + <9L dC l \ ... l[ .. T T ^ / _ T . 

- tt V 77" " 777 - 777 2 K « - F ^ u + M " 

+ (- J£ - J£F)u + (f + F£ - J£L) 



(8) 



Since the above equation is an identity equation with respect to u and u\ we obtain four 
differential equations. The four equations are readily solved to yield 

F = -coth{ l -(3i) , (9) 
J=-v1ete[sinh{±l3Z)]- 1 , (10) 
L = -ve-^sinhi 1 -^)]- 1 , (11) 



C = Tr 



±P£ + ln[sinh{^)] - v ] vcoth{^) 



(12) 



In the above, Grassmann variables v and are introduced as constants of integrations 
in order to make p a Gaussian type density matrix [8]. The normalized density matrix is 
obtained using these relations, 

p nor = l e Tr [-5^]£, e t[ sin ^(l^)] e «tt,-«t« e -(«tW)Mh(§/3€)](u+«) , (13) 

J V 

where the normalization constant H is obtained by the relation, 1 = Tr[p nor ], 

M = e Tr ^^Det[Acosh 1 -^} . (14) 

This normalized density operator corresponding to the above normalized density matrix is 
given by 

Pnor = Det[Acosh(\(3t)] 
The Helmholtz free energy is defined by 



e-P* . (15) 



(3 F =< lnp nor >+p <H> . (16) 



Substituting Eq.(15) into Eq.(16) and using Eq.(13) to calculate the expectation value, we 
obtain 

F = Nfi- k B T Tr[ln{\ + e - ^)] - Tr[ln2] . (17) 

The Fourier transformation of the above equation gives the familar result as follows, 

F = N/j, — k B T + e~ Kk )] ~ Nln2 , ( 18 ) 

where ^ k = e° k — /j, = — /j, and a is the spin index. The total energy is also readily 
obtained 

E = < H > 

= J2 DuDu ] h A Bi' ] A l l , BPuv\u=v 
A,B J 

= T,4n k , (19) 

fc,cr 

where Kab = £,ab + H and 

= tt^> • (20) 

These results on the non-interacting electron system clearly show that a finite temperature 
many-particle theory can be successfully formulated using the density matrix method based 
on the FSP approach. 

III. INTERACTING ELECTRON GAS SYSTEM 

Next, we investigate the interacting electron gas system of which the grand Hamiltonian 
is given by 

K = Y, Ub^b + \Y. VabVa^bMa , (21) 

A,B 1 A,B 

where Vab is any appropriate particle-particle interaction potential. The density matrix for 
an interacting system is not a Gaussian and generally quite complex. As a first approxi- 
mation, we choose a trial Gaussian density matrix with a parameter function, which will 



be chosen from minimization of the thermodynamic potential. Thus, the normalized trial 
density matrix is given by 

pnor ~ Det[Acoth{\(3Q)] 6 ' ' 1 ] 

where Q is the adjustible parameter matrix. The normalized density operator is written in 
the form 

pnor ~ Det[Acosh{\l3Q)} ' { ' 

in analogy with the free particle case. The thermodynamic potential is defined by 

Pfl = (3{F - fiN) 

= < lnp nor > +13 < H > -fifX < N > 

= < lnp nor > +13 < K > . (24) 



Substituting Eq.(23) into Eq.(24), we obtain 

pn = \$Y.Qaa - Y\H^osh{\pQ)\] AA 

1 A A 1 

— ~P ^ QabGba + ~P ^2 £abGba 

1 A,B 1 A,B 

+-(3 v ab[GaaG B b - GbaGab] , (25) 

H A,B 

where G = I — tanh(^PQ), and / is the identity matrix. The Fourier transformation of the 
above equation is 



pn= - V J dKtr[ln[4coshA(k)}} 
+ V j dKtr[A(k)tanhA(k)} 
l -VP J dKtrUk)G(k)} 



+ 

+ ^VPV(0) I <lKtr[G(k)] j (lK'tr[G(k') 



\vp f dKdK'V(k-k')tr\G{k)G{k')] , (26) 
8 J 
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where A(k) = ^(3Q(k), G(k) = a — tanh{\j3Q{k)) , tr is the trace for the spin indices only, 
and <T is the 2x2 identity matrix. Introducing a constraint parameter N, 

N = l -V [ dKtr[G(k)] , (27) 



2 

we rewrite the thermodynamic potential as follows, 

PS1 = a N - J dKtr[G{k)\ 
-V J dKtr[ln[4coshA(k)}} 
+V J dKtr[A(k)tanhA(k)} 
+ l -V(3 J dKtrUk)G{k)] 

--V/3 f dKdK'V(k - k')tr[G{k)G{k')\ , (28) 
8 J 

where a is the Lagrange's undetermined multiplier. Taking variations on j3Vt with respect 

— * 

to and A(k), we obtain 

« = - V -f?N , (29) 

PQ(k) = (fc(k) - a )a j dK'V{k - k')G{k') . (30) 
Eq.(30) is a self-consistent equation which determines Q(k). Since all functions in the above 

— * — * 

equation are 2x2 matrices, Q(k) and G(k) are also represented as linear combinations of 
the Pauli matrices. 

First, we consider the symmetry conserving paramagnetic case. In this case, the eigen- 

— * — * — * 

values of G(k) are same for both spins, and, thus, Q(k) and G(k) must be multiples of the 
2x2 identity matrix, 

Q(k) = q(k)a , (31) 



G{k) = g(k)a . (32) 
7 



Defining q(k) = e k — we can rewrite g(k) as follows 



g(k) = 1 - tanh[-/3(e k - //)] 



1 _|_ e P(e k -fj,) 

= 2n k . (33) 

Here, n k is the Fermi distribution function of the interacting system. The energy spectrum 
of the interacting system is obtained from Eq.(30) 

e k = e° k + (2S + l)V(0)n - J dK'V(k - k')n y , (34) 

where n = N/V is the number density The above equation is the familiar self-consistent 
Hartree-Fock equation which determines the energy spectrum of the interacting electron gas 
system Using the above results, we readily obtain the total energy and the Helmholtz 
free energy 

E=<H>= (2S + 1)V J dK(e° k + ^ k )n k , (35) 

and 

F = CI + Nfi 

= Nfi- Nk B Tln2 - k B T(2S + l)V J dKln[l + e - (ek -" )/fcflT ] , (36) 

where the self-energy is given by 

S fc = (25 + l)V(0)n - J dK'Vik - k')n k , . (37) 

Next, we consider the magnetic symmetry breaking solution. In order to represent the 
ferromagnetic state, we express Q{k) as follows 

Q(k) = (e k - n)a Q + 7 fc a 3 . (38) 

Assuming that ^ k is much smaller than e k and expanding to the first order of 7^, we obtain 



8 



tanh-8 



(e k -fi)a + j k a 3 



cr tanh^(3(e k - fi) + o^-^ k sedr?\^p{e k - fi) . (39) 



Using Eq.(30), Eq.(38) and (39), we obtain for the ferromagnetic state 



f3[(e k - fi)a + lk a 3 ] = - a)a - l - J dK'V(k - k') 



cr (l - tanh-(3(e k > - /i)) 



-o"3^7fc'sec/i 2 ^/3(e^ - jj) 



(40) 



From this identity equation, we obtain two self-consistent equations as follows, 

e ka = e + (2S + l)V(0)n - J dK'V(k - k')n k , a , (41) 

lk = \j dK'Vik - k') lk ,sech 2 l^[3(e k , - »)] . (42) 

Eq.(41) shows that the one-particle energy is now dependent on the spin index. Solving 
Eq.(42) self-consistently, one can obtain the amount of splitting between the up- and down- 
energy spectrum. Also Eq.(42) yields the celebrated Stoner condition as follows 

dn(e k 



1 = V(0) J dK 
= VF(0) 
« VN(0) . 



de k 



(43) 



Here, F(0), N(0) and V(0) represent the static Lindhard function, the density of state at 
the Fermi surface and the average exchange interaction respectively. The total energy and 
the thermodynamic potential are given by 







E = VJ2J dK[e° k + ^ k<T ]n ka , 

I _|_ g-/3(«fc<r-M) 



O = --In2-^VJ2 J dKln 



(44) 
(45) 



where 



S fc(T = ^V(0)n -±J dK'V(k - k')n k , a , 



n ka = 



I _|_ e /3(efc CT -/i) 

e ka and n ka can be determined self-consistently once j k is obtained through Eq.(42). 



(46) 



IV. CONCLUSION 



In this paper, we have presented a finite temperature many-particle theory based on 
the functional Schrodinger picture formalism. First, we have shown how the functional 
Schrodinger picture formalism can be applied to many-particle systems by solving the free- 
particle system exactly. Next, using the Gaussian density matrix obtained from the free- 
particle system, we have constructed a variational calculation formalism for the interacting 
electron gas system. The present variational formalism yields both the paramagnetic and 
the ferromagnetic solutions as in the case of the zero-temperature formalism. For the ferro- 
magnetic state, we have obtained a self-consistent equation which determines the splitting 
of energy between the two spin-split states. Also this equation leads to the Stoner criterion 
for the ferromagnetism. 

Here, we note that the above theoretical scheme is general and, thus, can be applied to 
any interacting condensed matter systems. Since the present theory is variational, and not 
perturbative, it is expected to become a useful tool for investigation of strongly correlated 
systems. In order to obtain results beyond the Hartree-Fock terms, it is necessary to carry 
out calculations beyond the Gaussian approximations. Variational methods employing non- 
Gaussian trial functionals and perturbational approach based on variational basis are under 
study and will be reported elsewhere. In conclusion, we believe that this application of the 
functional Schrodinger picture theory to finite temperature clearly demonstrates the versality 
and possible applicability of the theory to wide range of condensed matter problems. 
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